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A new analytical model is developed for the prediction of noise from serrated trailing- 
edges. The model generalizes Amiet’s trailing-edge noise theory to sawtooth trailing- 
edges, resulting in an inhomogeneous partial differential equation. The equation is then 
solved by means of a Fourier expansion technique combined with an iterative procedure. 
The solution is validated through comparison with finite element method for a variety 
of serrations at different Mach numbers. Results obtained using the new model predict 
noise reduction of up to 10 dB at 90° above the trailing-edge, which is more realistic 
than predictions based on Howe’s model and also more consistent with experimental 
observations. A thorough analytical and numerical analysis of the physical mechanism 
is carried out and suggests that the noise reduction due to serration originates primar¬ 
ily from interference effects near the trailing-edge. A closer inspection of the proposed 
mathematical model has led to the development of two criteria for the effectiveness of 
the trailing-edge serrations, consistent but more general than those proposed by Howe. 
While experimental investigations often focus on noise reduction at ninety degrees above 
the trailing-edge, the new analytical model shows that the destructive interference scat¬ 
tering effects due to the serrations cause significant noise reduction at large polar angles, 
near the leading edge. It has also been observed that serrations can significantly change 
the directivity characteristics of the aerofoil at high frequencies and even lead to noise 
increase at high Mach numbers. 


1. Introduction 

The past few decades has seen a rapid growth of air traffic, while the public’s attention 
to aircraft noise and its health consequences has also been continuously increasing. This 


has led to more stringent regulations for aircraft noise (Casalino et al. 2008). With 


regard to the impact of aircraft on community noise, the take-off and landing process are 
of main concern. Among the different mechanisms present during the landing process, 
airframe noise is believed to be the dominant component. It is widely accepted that the 
broadband noise, induced by the interaction of boundary layer with the aerofoil trailing- 
edge, known as the turbulent boundary layer trailing-edge noise, plays a significant role 
in the overall airframe noise. Turbulent boundary layer trailing-edge noise also dominates 


the noise produced by wind turbines (Oerlemans et al. 2007). Unless explicitly stated. 


the turbulent boundary layer trailing-edge noise will be referred to as trailing-edge noise 
in the rest of this paper. 
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Figure 1. The schematic of a fiat plate with and trailing-edge serrations. 


When a turbulent boundary layer convects past the trailing-edge, unsteady pressure 
with a wavenumber in the hydrodynamic range is scattered into sound (Chase 1975). 
Both experiments and theory reveal that the radiated sound power varies with the flow 
velocity to the power of 5, which is more efficient, at low Mach numbers, compared to 
the power of 8 valid for free stream flows ( LighthilT||1952 Williams fc Hall||1970 ). 

Different models have been put forward for predicting and understanding trailing-edge 
noise. In 1976, Amiet (19766 19781 proposed a semi-analytical model in which the aerofoil 
is modelled as a flat plate. The model followed Schwarzschild’s technique to obtain the 
scattered pressure on the surface of the flat plate and the far-held sound was obtained 
using the surface pressure integral based on the theories of Kirchoff and Curie. The 
model established an analytical relationship between the far-held sound spectral density 
and the wavenumber spectral density of the wall surface pressure under the turbulent 
boundary layer. Amiet’s model agrees well with experimental observations, especially at 
high frequencies ( Roger fc Moreau|2005 1. It is worth noting that Amiet’s model assumed 
that the leading edge of the plate is inhnitely far away from the trailing-edge and thus 
has no effects on the scattered pressure, which might not be accurate at low frequencies. 
In 2005, in order to investigate the leading-edge back-scattering effects, [Roger fc Moreau 


(2005) extended Amiet’s model by incorporating the back-scattered pressure from the 
leading edge and found that when the Helmholtz number kc > 1, the back-scattering 
can be safely ignored and only at very low frequencies does the back-scattering alter the 
far-held sound. 

As trailing-edge noise dominates the sound generation at low Mach numbers, different 
noise reduction techniques have been investigated. Howe proposed a theoretical model 
to predict the sound generated by a semi-inhnite plate with serrated trailing-edge of 
sinusoidal and sawtooth prohles ( jHowe|[l991a|6 ), see hgure[^ Howe’s models show that 
sawtooth serrations are more effective in reducing the trailing-edge noise than the si¬ 
nusoidal ones and that the use of sharp sawtooth serrations, i.e. 2h/X > 8, can lead 
to signihcant reduction of trailing-edge noise. More recently, Azarpeyvand et al. (2013) 
carried out an analytical investigation of trailing-edge noise reduction using novel serra¬ 
tions, namely, sawtooth, sinusoidal, slitted, slitted-sawtooth and sawtooth-sinusoidal. It 
has been found that the noise reduction is a sensitive function of the complexity of the 
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serration geometry and significant noise reduction can be achieved by applying complex 
periodic serrations to the trailing-edge. It has also been shown that the slitted-sawtooth 
serration is the most effective design amongst the aforementioned serration geometries. 
An experimental investigation on trailing-edge serrations was performed by [Dassen 


et al. (19961. Both aerofoils and flat plates of different shapes were tested in a wind 


tunnel. A maximum of 10 dB noise reduction for the flat plate and 8 dB reduction for 
the aerofoils were reported, both of which occurred mainly at low frequencies. Later, 


Parchen et al. { 19991 conducted an experimental investigation of the aeroacoustic effects 


of trailing-edge serrations on wind turbine blades at both full and wind-tunnel scales. An 
average sound reduction slightly below that reported by Dassen et al. was observed. Most 


recently, Oerlemans et al. (20091 examined and compared the noise generated by stan¬ 


dard, serrated and aeroacoustically-optimized aerofoils, but only 2-3 dB noise reduction 
was achieved for the aerofoil with serrated trailing-edge at low frequencies. Both |Parchen| 
et al. (1999) and Oerlemans et al. (2009) reported a noise increase at high frequencies. 


Gruber (20121 recently conducted an extensive experimental investigation on the aeroa¬ 


coustic performance of aerofoils fitted with different sawtooth and novel serrations. The 
acoustic measurements were performed to give the sound power level (SWL) integrated 
in the mid-span plane. An average of 3-5 dB reduction was achieved using sharp sawtooth 
serrations, and a noise increase of up to 5 dB at higher frequencies was also reported. 
It was explained that the significant reduction of phase speed near the sawtooth edges, 
together with a slight reduction of the coherence of pressure measured along the edge 
is responsible for the sound reduction observed in experiments. All the experimental 
studies, however, indicate that Howe’s model significantly overpredicts the sound re¬ 
duction capability of trailing-edge serrations. This might be caused by the assumptions 
and approximations used in Howe’s derivation, which will be discussed in detail in the 
subsequent sections. 

Although different serrations have been used in many applications such as wind tur¬ 


bines and jet nozzles (Oerlemans et al. 2009 Callender et al. 2005 Yan et al. 2007), 


the physical mechanism of the noise reduction remains poorly understood. Howe’s model 
gave a first insight into the physics involved, but the large deviation from experiments 
suggests a new and more accurate theory is needed. 

The main objectives of this paper are to develop a new theory to predict the sound 
generated from a serrated trailing-edge more accurately, and to improve our understand¬ 
ing of the sound reduction mechanism due to the presence of trailing-edge serrations. 
The paper is structured as follows: the new theoretical model for sound radiation from 
serrated trailing-edges is presented in section Section provides validation results 
against finite element solution. A parametric study will also be presented and the effects 
of serrations on trailing-edge noise directivity will also be discussed. A comparison be¬ 
tween the new model and Howe’s model is made in section The physical mechanism 
of sound reduction using trailing-edge serrations is discussed in section and noise re¬ 
duction criteria developed based on the results in section will be listed and discussed. 
A brief conclusion is given in the last section. 


2. Analytical formulation 

As for Amiet’s model, the analytical model developed here is based on Schwarzschild’s 
technique for wave equation with a discontinuous boundary condition. It is therefore 
useful to begin with describing this technique. The Schwarzschild method ( Amiet|19766 
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Roger fc Moreau[2005 ) states that if a function f{x,y) satisfies 

d^f , 52 / 


then, for a; < 0, 


|(x,0) = 0, 

, fix, 0 ) = g{x), 


a; < 0 
a: ^ 0, 


( 2 . 1 ) 


1 /■°° 


( 2 . 2 ) 


As shown by Amiet ( 1976a|), the above method can be used to obtain the scattered 


pressure field over the surface of the aerofoil. 


2.1. The mathematical model 

Consider an aerofoil with trailing-edge serrations, modeled as a flat plate as shown in 
figure with an infinitesimal thickness and an averaged chord length c and spanwise 
length d. Let x', y' and z' denote the streamwise, spanwise and normal to the plate 
coordinates, respectively. The observer point is located at (xi,X 2 , X 3 ). The profile function 
H{y') is used to describe the serrated edges. The origin of the coordinates is chosen in 
such a way that H[y') is an oscillatory function of zero mean and that H{y') = 0 in the 
absence of serrations. Figure [T] shows a sawtooth serration with a root-to-tip length of 
2h and a periodic wavelength of A. 

When the sound wavelength is equal to or shorter than the chord length c, the flat 
plate can be treated as a semi-infinite plate without a leading-edge ( Amiet|[T976fe I. Fur¬ 
thermore, the plate can be considered infinite in the spanwise direction provided it has a 
relatively large aspect ratio (typically d/c > 3) ( Amiet|1978 Roger &: Carazo|2010 \ . The 
turbulence inside the boundary layer is assumed to be frozen, i.e. it remains statistically 
the same before and after passing over the trailing-edge. 

After implementing a spatial and time Fourier transformation, the hypothetical surface 
pressure beneath the turbulent boundary layer that would exist when the flat plate is 
absent can be expressed as an integral of different wall pressure gust components. The 
incoming wall pressure gust of frequency w convected at a speed Uc, as illustrated in 
figure takes the form of 


-k-2v') 


(2.3) 


where Pi is the magnitude of the incident wall pressure gust and ki and k 2 denote the 
wavenumbers in the chordwise and spanwise directions, respectively. 


The sound sources due to the presence of solid boundaries (Curie 19551 can be modelled 


as dipoles, in addition to the quadrupoles in free field (Lighthill 19521. As explained in 
Amiet’s paper ( Amiet|1976b I, the incident pressure produces a scattered field originating 
from the trailing-edge, due to the change in boundary condition at the wall. The scattered 
held induces a pressure jump that cancels the incident pressure jump at the trailing- 
edge and in the wake after the plate (Kutta condition). Thus, the total pressure can be 
decomposed into two parts, namely pt = Pi p. The incident wall pressure is given by 
2.3 and the scattered pressure field, p, must satisfy the following conditions at z' = 0, 


dp 


= 0 , 




. = _p.^-d^t-k^T-k^v') ^ ^ 


(2.4) 
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It is worth pointing out that the incident pressure pi and the scattered pressure p on 
the surface of the plate are, as in Amiet’s paper ( Amiet|[l976& ), identical to the incident 
pressure jump and the scattered pressure jump, respectively, as mentioned in the pre¬ 
ceding paragraph. This is merely a mathematical simplification and can be shown to be 
valid through a more rigorous analysis. 

In the plate-fixed frame {x', y', z'}, the air flow has a uniform speed U in the streamwise 
direction outside the boundary layer and the wave equation governing the scattered 
pressure field p is 


I 


d 


d 




(2.5) 


^ dt dx' 

where cq denotes the speed of sound. With the assumption of harmonic perturbation 
p = P(x', y', 2 :')e'‘“‘, the above equation reduces to 




( 2 . 6 ) 


where fc = w/cq, = I — Mg and Mq = U /cq. 

In order to make the boundary conditions in (2.4) independent of y', the coordinate 
transformation (Roger et al. 2013) x = x' — H{y'), y = y', z = z' is used and leads to 
the following differential equation: 

B'^P BP 

+ ^ = “• 

(2.7) 

where H'{y) and H"{y) denote the first and second derivatives of H{y). The boundary 
conditions now read 


[P^ + H'\y)) 


V(x,y,0) = ^ ^ ^ Q 

BP{x,y,Q)/dz = 0, x < 0. 


( 2 . 8 ) 


Since the coefficients in (2.7) are y-dependent, the standard “separation of variables” 
technique cannot be applied to solve this equation. We therefore turn to using a Fourier 
expansion technique in the following derivation. 

2.2. Fourier expansion 

As the scattering problem is periodic in the spanwise direction, one can expand the 
scattered pressure field using Fourier series as 


P(x, y, z) = ^ P„(x, z)e*'"""^, 


(2.9) 


where k 2 n = ^2 + 2n7r/A. 

Substituting the above expression into the transformed wave equation, shown in (2.7), 
yields 


(^3^ + H'\y)) 


5x2 


5y2 


52 02 


(2iMoA: - H"(y)) 


^P„(x,z)e‘ 


ik-inV 


= 0 . 


( 2 . 10 ) 


Multiplying (2.10) by then integrating it over y from —A/2 to A/2, it can be 
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Figure 2. The schematic of sawtooth serrations 


readily shown that 


~ ^ 2 n') ^ Pn'{x, z) + 


dx"^ 

1 

^ " -•/ - n——oo ' 

( 2 . 11 ) 

Note that when both H'{y) and H”{y) are constant within an entire sawtooth wave- 


OO 


- (-ff"(2/) + 2ife„-ff'(y))^|p„(x,0)e*[2("-')-/A]i/dy = 0. 


length, the summation over different modes in (2.11) can be dropped and we obtain a 


fully decoupled differential equation for mode n'. However, this means that the flat plate 
has a straight or swept trailing-edge. For serrations of an arbitrary profile, both H'{y) 


and H"{y) generally depend on y. Thus 2.11 becomes a coupled differential equation, i.e. 


more than one mode appears in each differential equation. The physical interpretation 
of this mode coupling will be discussed later. 

In this paper, we only focus on the sawtooth serration, which has been shown to be 


effective in reducing the trailing-edge noise Howe (1991&). The method can however also 


be used for other serrations. Consider a sawtooth centred around the coordinate origin, 
and let (xoj eq)) (Xi) Ci) (x 2 , ^ 2 ) be the Cartesian coordinates of the tip and roots of 
the sawtooth, as shown in figure The serration profile function H (y) can therefore be 
defined as 


H{y) 


0 - 0 ( 2 / - Xo - mX) -f eo, 
0 - 1 ( 2 / - Xi - "1^) + eu 


Xo + wA < 2 / ^ Xi -I- mA 
Xi -I- toA < 2 / ^ X 2 + rnX, 


( 2 . 12 ) 


where cTj = (e^+i — e^)/(xt+i ~ Xi)i J = 0; 1 and m = 0, ±1, ±2, ±3 • • •. Let a, defined as 
(7 = \aj\ = 4h/A, denote the sharpness of the sawtooth serrations. For a sawtooth profile, 
as mentioned above, H'{y) is not continuous and H"{y) is thus singular at the joint- 
points. We use the conventional generalized function S(x) to describe the singularities at 
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these points, i.e. 


H'(.y) = 


CTo, Xo + wA < y < xi + wA 
(Ti, xi + ^7^A < y < X2 + wA 


H''{y)= ^ {-l)'^+^2a5{x-m\/2). 


(2.13) 


m——oo 


Substituting the serration profile function and its derivatives, ( |2.12[ ) and ( |2.13[ ), into 
the wave equation, (2.11), and making use of the fact that f(x)S(x — T)dx = fir), 
we obtain 


^2 ^“2 Q 

(/3^ + cr^) ^ 2ikMo— + (/c^ - kl„,) )■ Pn'{x,z) 


4(7 

T ^ 


i—n'—odd 


1 - 


fc 2 A + 2n7rA dPn{x,z) 


(n — n')'K 


(2.14) 


dx 


To make the above equation more compact, we write the set of differential equations 
obtained above into a matrix form. Let a linear operator 




and a vector of functions 

P — (’ ’ ’ P—n' (x^ zf P—n' + l {x, Z^ j • • • Pji' — I (x, z), Pji' (x, z), ‘ * * ) 
then the coupled equations can be written as 


I)P-AP=B 


dP 


(2.15) 

(2.16) 

(2.17) 


where the symbol T in (2.16) represents the transpose of a matrix. Matrices A and B 
denote the coefficient matrices of P and dP/dx, respectively. Ami and Bmi representing 
the entry corresponding to mode m in row and I in column of matrices A and B are 

e 4(7 TO + Z + k2X/TT 

I-;-, TO — / IS odd 

l-m (2.18) 

m — I is even. 


^ml — (^2m ^ )dmh Pml — ^ X 

lo. 


where Smi represents the Kronecker delta. It is worth noting that the index of matrices 
A and B are from —to to to and —I to I rather than from 1 to 2m + 1 and 1 to 2Z + 1 
given the fact that the mode number are symmetric with respect to 0. 

Substituting the profile geometry, (2.12), into the boundary conditions, (2.8), and 
performing the same Fourier expansions, yields 

Pn{x, 0) = x> 0 

X ^ 0, 




(2,19) 


where a„ are defined as 


f .\/2 


a„ = - 


-A/2 




( 2 . 20 ) 


In (2.17), as B is not a diagonal matrix, the term B{dP/dx) contains coupling terms, 
in the sense that different modes, for example, appear in the governing equation 
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of Pm- This means that every mode is interacting with the other modes and cannot 
be solved independently. From the expression of B in (2.18), it can be seen that the 


strength of coupling is propotional to a/X. This indicates that shaper serrations have 
stronger coupling between different modes. Also note that A is a diagonal matrix, so if 
6 « 0, i.e. the serrations are very wide, then there is no coupling effect and one can solve 
each mode independently. 

At very low frequencies, the contribution of higher modes is expected to become grad¬ 
ually negligible compared to mode 0. Thus, it can be reasonably assumed that in the 
case of the governing equation of Pq (mode 0), the coupling with higher modes is weak 
and one can solve Pq individually. The coupling effect becomes more pronounced at high 
frequencies and for sharp serrations. To solve these coupled equations at relatively high 
frequencies, one can use an iterative procedure, to be introduced in the next section. 

2.3. The iterative solution procedure 
2.3.1. Scattered surface pressure 


To obtain the scattered surface pressure, equation (2.17) together with the boundary 


conditions in (2.19) need to be solved. By analogy with the solution of a system of linear 
algebraic equation ( Siili fc Mayers|2003 ), which can be obtained iteratively, we shall solve 
our system of PDFs in an iterative manner. 


Substituting a known initial value into the coupling term in (2.17), one can write 

dp(o) 


DP - AP = B- 


dx 


( 2 . 21 ) 


Solving (2.21) yields a new set of solutions Replacing P^^'> in (2.21) with P^^\ we 

dpw 


obtain a new wave equation. 


DP - AP = B- 


dx 


( 2 . 22 ) 


Solving (2.22) gives a new set of solutions P^^\ This process is repeated to obtain a 
solution sequence, P^^\ P^^\ P*^^\ P^^'> • • •. If the sequence is convergent, then its limit 


satisfies ( |2.17[ ). 

The initial value P^°^ used to start the first iteration can be obtained, as mentioned in 
the previous section, by ignoring all the coupling terms, i.e. with 6 = 0, and by solving 


each equation individually via the standard Schwarzschild technique, as shown in (2.2). 
The solution to each equation in the decoupled system of equations 


DP - AP = 0 


(2.23) 


can be found as follows. 

After use is made of the transformation of Pn' = individual equa¬ 


tions in (2.23) reduce to. 




^ +O-)}P„'-0, 


dx"^ dz"^ 


where 


Kr, 


= yfc2(l + a2)-fc|„,(/32 + a2)/(/32 + a^). 


(2.24) 


(2.25) 


Making use oi X = x, Z = + a'^z, one can show that (2.24) reduces to a standard 

Schwarzschild problem and the solution can be found using the Schwarzschild integral 


described in (2.2), as 


pW = P,e''=i"a„-((1 - im-Tn'x) - 1), 


( 2 . 26 ) 
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where is the element of the vector corresponding to the n'-th mode, and 


l^n' — Kn' + fci + 


kMn 


/32 


r e'* 

Etx) = / , dt. 

Jo V^t 


(2.27) 


The initial solutions obtained by ignoring all the coupling terms denote the decoupled 
part of the exact solution of each mode, which implies that the n-th mode excitation 
{x > 0) produces only an n-th mode response (a; < 0). The iteration procedure will add 
a coupled part to the solution of each mode. The coupled part implies that an n-th mode 
input {x > 0) will also produce some m-th mode responses {x < 0), where m 7 ^ n. It 
can be expected that those coupling contributions from closer modes are stronger than 


remote ones. For sawtooth serrations, as will be shown in (2.30), the coupling strength 
decays quadratically with respect to the difference between their “mode number”, i.e. 
\m — l\. 


By substituting into the coupling terms on the right hand side of (2.17), one 


obtains some inhomogeneous equations that can no longer be solved using the standard 
Schwarzschild technique. However, if one can transform these equations into homog eneou s 
ones, then Schwarzschild’s method can again be applied. Note that satisfies (2.23), 
hence, for cc 7 ^ 0 , where is first-order continuously differentiable, the following 

equation holds: 


5p(o) 9p(o) 

D—-/\- 


dx 


dx 


= 0 . 


Making use of (2.28), (2.21) can be equivalently written as 


D P 




where v is a coefficient matrix whose entries are 

. - 4/1 

Pral I 

— 


^2 _ 1.2 
^2m ^21 


= z 7r2(m — ly 


spm\ 


m — I = odd 


(2.28) 


(2.29) 


(2.30) 


0 , 


.-I = 


It is worth pointing out that (2.29) only holds when a; S R and a: 7 ^ 0, and in order to 


apply the Schwarzschild technique, it must be valid over the whole domain. However, since 
the singularity of dP^^^dx only exists at x = 0 , similar to the differentiation of H{y), 
we may again make use of the generalized function to account for this singularity. Let 
dP^^l/dx denote the generalized differentiation, which allows the presence of generalized 
functions at singular point x = 0 but equals to dP^^'>/dx elsewhere, then equation 


5p(o) 5p(o) 

D— -/A—— = 0 


dx 


dx 


(2.31) 


needs to hold over x S R. The Schwarzschild technique suggests that if ( |2.31[ ) does hold, 
then the routine application of the steps described from (2.24) to ( |2.26[ ) shall recover the 
value of /dx for x < 0. Thus, one can verify that the intended dP^^^dx can indeed 
be found as 


dP^°'> 

^ n' 

dx 


dP^P 


(a:, 0 ) = (x, 0 )-f P»an'(l - i){-y/JI/f )Vp7rx6{x), 


(2.32) 
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where /dx denotes the element of dP^^pdx, corresponding to the n'-th mode and 


1 

S{x)dx=-. 


Now, the first iterated solution can be obtained by solving equation 

dpio) \ 


D P 


dx J 


A P 


dpio) \ 


(2.33) 


(2.34) 


via the ste ps de scribed from (2.24) to (2.26). 

Solving (2.34) gives the values of Continuing this iteration process gives P^‘^\ 
P^^'> ■ ■ ■. The exact solutions P can also be expressed as 

P{x, 0) =N{x) + (x) + (x) + C^^'> (x) + • • • , (2.35) 

where N is the non-coupled part and the coupled parts are denoted by = P^®^ — 
p(*-i) (j = 1, 2, 3 • • • ). The entries of N and corresponding to mode n' are given 
by 


Nn' (x) = P^e 


Cib(x) = Pe''=^"(l-i) ^ 


ikix ^ 


((1 - i)E(-/x„/x) - 1), 


m^m 


m— — oo 


^ifci(^E(-/x„/x) - E(-/x,„x)^ 

V^(' 


(2.36) 


(2.37) 




The elements of the second-order function are provided in appendix A. Since Vmi oc h 

(i) 

it can be readily shown that C^, oc h* for i = 1, 2, 3 • • •. This means that the solution 


presented in (2.35) is a perturbation (Taylor-expansion) series with respect to half of the 


root-to-tip length h. Therefore a smaller value of h compared to the sound wavelength 
yields faster convergence. It is worth noting that the function becomes more and 
more complex as i increases. However, if vanishes sufficiently quickly, higher orders 
can be neglected without causing significant errors. This appears to be the case for the 


frequencies relevant to trailing-edge noise, see Sec. 3.2 Substituting (2.36), (2.37) and 


(A2) into (2.35), a second-order approximation of the exact solutions is obtained. 


The scattered surface pressure is obtained by summing P„/ (x, 0) over all different 
modes and transforming back to the physical coordinate system, namely 

OO 

P(x',2/',0)= ^ P„.(x'-P(y'),0)e‘'=-'^'. (2.38) 

n'——oo 

Here, P„' is the solution obtained from the iteration procedure mentioned above, 

OO 

P(x',j/',0)= ^ + ---){x'-H{y'),0)P'^^^'y', (2.39) 


n — — OO 


where Nn' and are defined in (2.36) and (2.37), and can be found in appendix 


A. Note that the terms in the second parenthesis are the arguments for the Nn' and 
C^) {i = 1, 2, 3, • • •) functions. It is worth pointing out that in the l imiting case w hen 
H{y') = 0, Cb) vanishes, and (2.39) reduces to the result obtained by Amiet (19766) for 
a straight edge. 


As shown in (2.39), the scattered pressure field can now be expressed in terms of an 


infinite series. By inspection of (2.39), one can show that at sufficiently low frequencies. 
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i.e. kih < 7r^/4, the infinite series is absolutely convergent. At higher frequencies, the 
series still appears to be convergent, but to obtain satisfactory approximation, a higher 
truncation number and higher-order iterations may be required. The convergence of the 
series will be disscussed in the following sections by comparing the far-field sound prediced 
using different order approximations. 


2.3.2. Far-field sound pressure 

As illustrated in figure[^ the observer point is located at (xi, X 2 , xa) and the flat plate 
has an averaged chord length c and span length d. The far field sound can be found using 
the surface pressure integral as mentioned in Amiet’s model (Lamb |1932 Curie 1955 
|Amiet|[T975t : 


p/(j 


— ICJX 3 


AP(x',?/')e'*''^dx' dy', 

where AP = P denotes the pressure jump, and Sq = xf -I- ^^(x 2 -|- x§) and 

Mq{xi-x')-So xix'+ X 2 y' 

^2 + ■ 


(2.40) 


(2.41) 


By substituting the solution obtained in (2.39) into (2.40), the far-field sound pressure 
can be expressed as 


, ,, ,^-ia;x3C^^ ^sin((A-f l/2)A(fc2-fcx2/5'o)) 

(ii^ j ^ si„(A/2fe-fa,/5„)) 

Here, 2N -\- 1 represents the number of sawteeth on the edge and the far-field sound 
gust-response function C is defined as 


-ik/l}'^{Moxt-So)Jk/p'^{Mo-xt/So)h 


£{uj, /c 2 ) — (1-i)— ( ^ (0„' -I- -I- -I-) j e'*' 

Vn'=-oo / 


with (only the first two terms are given, see more results in appendix A) 

0 n' — ^n'Qn'n'^ 

00 

0^/ — ^ ^ ^n'ndk\arri{f^n'n' Qn'm) rn'rny/ t^n' Sn'm'j- 


(2.43) 


(2.44) 


m——oo 


The functions Qnm and Snm in the above equations are given by 
1 


j=o 


Qnm = ^ (e‘”"’^^+ir(c-b ej+i;^™,^^) - e*''""^T(c-b ^a)) 

- ^nj ^ f-^A 


g njiXj C + j)/ + Cj-i-i 5/i-m,,) \^{c-\- IXrn-) l^Bnj 

l^Bnj 

(2.45) 

1 , , 

Snm = V ^-( . (e‘”"’^"+"E(?7Am(c + Cj+i)) - e'''""£( 77 ^™(c + ej)))) 

^ IKnj V .yVAm 

- -^e--(x.-(c+.H/-H(E(,7B™^.(c+e,+0) -E(r7B„^,(c + e,^ 
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where the function F is defined by 


r(a;; v) = e‘“"'^E(^a;) 



- iz)x) + 


1-i 




(2.46) 


fij^ = ki+ k/p'^{Mo - xi/So), 

f-^Bnj — (^2n kX2lS q)/ fJj , 

Knj = k2n “ kx 2 /So + k/{M q - xi/So)aj, (2.47) 

VAm. = K^ + kMo/dd"^ + cr^) - k/0^{MQ - xi/Sq), 

VBmj = Kjn + kMo/{l3'^ + a^) + (k 2 „ - kx 2 /So)/aj. 


Note that (2.42) is the far-held sound induced by the scattered pressure only. When the 
incident pressure is also incorporated, as pointed out by Amiet (1978), the number 1 
appearing in the parenthesis of function F should be omitted, i.e. the third term on the 
right hand side of (2.46) should be replaced by —e'‘‘"“/(l — i). 


2.3.3. Statistical formulation 

The hypothetical surface pressure of frequency oj beneath a turbulent boundary layer 
on the plate surface that would exist when the semi-inhnite plate is absent can be ex¬ 
pressed as a Fourier integral. 


X, y')= Jl fci, dfci dfc2. (2.48) 

Generally, for a given frequency w, ki can have different values ( Amiet|19766 ). However, 
experiments (Willmarth 1959) have shown that Pi{u!,ki,k 2 ) peaks in the vicinity of 
ki = ujfUc, where the convection velocity Uc is only a weak function of ui. Hence (2.48) 
reduces to 


/ OO 

P,ioj,k2)P^'^^^'+^^y''> dk2. 

-OO 


(2.49) 


As shown in the preceding section, a wall pressure gust of 


will induce a far-field sound pressure 

f -iuJX 3 c \ ^ sm{{N+l/2)X{k2-kx2/So)) 

sm{X/2{k2-kx2/So)) 2) *(w, 2). (■ ) 

Thus, the wall pressure defined by ( |2.49[ ) will induce a far-field sound pressure of 

, , f-iujX 3 c\ sin {{N+1/2)\{k 2 -kx 2 /So)) 

J.J sin(X/ 2 (M-k../S,)) 

(2.51) 

The PSD of the far-field sound is given by 

Spp{x,oj) = ^li^ (||;(p/(x,w)p}(x,w))) , (2.52) 


where the asterisk denotes complex conjugate, and 2T is the time length used to obtain 
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Pf{x,u}) by performing Fourier transformation. Substituting ( |2.51[ ) into (2.521 yields 

/ UIX 3 C 


S'pp(x,a;) = 




( sin ((TV + l/2)A(fc2 - fcxa/^p)) 
, i sin (A/ 2 (A ;2 - fca:2/5'o)) 


|/:rn(u;,fe)dfc2, 

__ (2.53) 

where n(w,A: 2 ) is the wavenumber spectral density (Amiet 19751 of the hypothetical 
wall pressure beneath the turbulent boundary layer on the plate surface. For very wide 
serrations, i.e. /i « 0, (2.53) reduces to Amiet’s model (Amiet 19766). Equation ( |2.53 ) 
can be simplified by assuming a very large span, i.e. the number of serrations (2N + 1) 
is sufficiently large. Using the following equation. 


lim A^ 
N—>-oo 


sin^ {{N + l/2)A(fc2 - kx2/So)) 
sin^ (A/2(fc2 - kx2/So)) 


27rd S{k 2 — kx 2 /So2rmr/X), (2.54) 


where S{x) is the conventional generalized function defined in Sec.[^ one can show that 
the PSD of the far-held sound in the plane y' = 0 is given by 

S'pp(x, w) = ( -— |£(a;, 2m7r/A)|^ n(a;, 2m7r/A). (2.55) 

' U Q / m= — f-v-i 


Equation 2.55 is the fundamental result of this paper and it is interesting to note that 


the inhnite series in (2.55) appears similar to that in Howe’s model shown in (4.2). For 


example, both results show that the PSD of far-held sound is related to the wavenumber 
spectral density of the surface pressure through n(a;, 2m7r/A), therefore a skewed wall 
pressure gust with k 2 = 2 m'KfX plays an important role in sound generation. 

2.4. Discussion on the effects of serration geometry 


The complicated formulation of the far-held noise (2.55), and the response function (2.43), 


make it very difficult to assess the effectiveness of serrations without numerical evaluation 
of the equations. This section attempts to derive two simple conditions for serrations to 
obtain effective noise reduction. In order to achieve signihcant sound reduction, we wish to 
minimize (2.55). Since \C{ui, 2m7r/A)p is very complex, we will perform an order analysis 
hrst. 

Careful examination of (2.43) shows that |£p is proportional to For illus¬ 


tration purposes, we assume that the observer is at 90° above the trailing-edge in the 
mid-span plane, i.e. xi = 0 and X 2 = 0, and that the Mach number is low, e.g. Mo < 0.2. 
Then Kmj ~ 0 when m satishes k 2 + 2m7r/A « 0 and thus the value of |£(a;, —2m7r/A)p 
is dominated by mode to, i.e. 

2 


iQmm + ’ 


(2.56) 


A2c2 

= 46,/(7r2(n — to)^) when n — to is odd, equation (2.44) 


\L(u}, —2TO7r/A)| 

Furthermore, noting that 
suggests that 0 m may roughly be approximated by only summing over modes to — 
1 and TO -|- 1, since higher orders to ± {2j 1) with j ^ 1 are at least one order of 

magnitude smaller due to the quadratic term in the denominator of Vnm- Usin g this 
approximation, ©m varies linearly with 0^-1 and am+i- It can be shown fr om ( 2.58) 
that |am-i| and |am-i-i| are of the order of \am\ and |©m | = 0(|am|)- From (2.56), we 
hence have |£(w, —2TO7r/A)p = 0(|amP), and 


|£(a;,2TO7r/A)pn(w,2TO7r/A) = O ( Y 2mTT jX) 


(2.57) 


m— — oo 
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We are now in a position to discuss the conditions for minimizing (2.57). From the 
definition of am in (2.20), it can be shown that 


t/2 


-sinc(fci/i — m'K/2) 


g—im7r/2 


sinc(fci/i + TO7r/2). 


(2.58) 


2 " '' 2 
Thus, \am\ is maximum when m « ir'o, where vq = 2kih/TT. To minimize the right hand 


side in (2.57), we therefore require that n(w, 2m7r/A) <C n(w, 0) when m approaches ±izq. 

(2.59) 


Assuming frozen turbulence, n(a;,2m7r/A) is given by 

2 TTm\ 1 


n UJ 


A / 27r 


Sqq(w,y')e 


dy'. 


When 2'KVQlyi/\ = kilyia ^ 1, the integrand in (2.59) for m close to (or larger than) 


t'o oscillates rapidly within the length scale ly> of Sq^ui, y'), which corresponds to the 
spanwise correlation length given by 

ly,{uj) = 


SqqipJ, 0 ) J —c 


Sqq{uj,y)dy. 


(2.60) 


The integral in (2.59) therefore evaluates to a small value compared to n(a;, 0 ). Thus, 
a condition for noise reduction is that kihe ^ 1, where we have defined an effective 
root-to-tip length 2 he = aly' that describes the correlated serration height. 

As 2kihg = kialy'ipj), the decay rate of spanwise correlation length ly'{uj) with respect 


to frequency is critical. If one makes use of Corcos’s correlation length model (Corcos 
1964), ly'{u!) « 2.1Uc/ui, kihe will reduce to a constant, in this case 2.1a, a sole function 


of the serration sharpness factor and independent of frequency. This is consistent with 
the findings of Howe (1991&). However, if the decay rate were faster than that given by 


Corcos, no sound reduction or even some sound increase would occur at high frequencies. 
An accurate description of the characteristics of the surface pressure fluctuation beneath a 
boundary layer is therefore critical for the model to accurately predict the sound reduction 
at high frequencies. 

Note that the condition kihg ^ 1 is only a necessary condition, because when kih ^ 0 


there is no noise reduction. This can be seen from equation (2.58), since am = 0 when 
kih —)• 0 except when m = 0, so the right hand side in equation (2.57) reduces to n(a;,0) 


which corresponds to the straight edge case. From (2.58), for a given integer m away 
from izq, \am\ tends to zero when kih^ 1. This provides another necessary condition for 
noise reduction. Physically, the height of the serrations must be sufficient for them to be 
seen by the incoming hydrodynamic waves. 

We have thus obtained two necessary conditions for noise reduction, kih^ ^ 1 and 
kih ^ 1, that are consistent with those proposed by Howe. These conditions will be 
further investigated in section 


3. Results 


3.1. Model Validation 


3.1.1. FEM implementation 

For the coupled differential equations mentioned in the last section, the solutions are 
obtained by performing an iterative-solving procedure. In this section, we shall investigate 
the validity of the proposed iterative solution using the Finite Element Method (FEM). 
Instead of solving the far-held sound directly, costing a signihcant amount of computer 
memory, a feasible alternative is to calculate the near held using FEM and obtain the far- 


held solution by performing a surface integral, as adopted in analytical models, see (2.40). 
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Figure 3. Illustration of the FEM computing domain and the boundary conditions 


In order to make a direct comparison between the computational and analytical results, 


the wave equation together with the boundary conditions given in (2.5) and (2.4), respec¬ 


tively, will be solved. The governing equation and boundary conditions in the frequency 
domain can be written as, 


^Q2p Q2p Q2p 


+ 


+ 


dx''^ dy'^ dz'"^ 


dP 

2ikMQ--—- -\- P = 0 , 

ox' 


^{x',y',0)=0, 

, P{x\y\0) = -P,e 


i{kix'+k2y') 


x' < H{y') 
x' > H{y'). 


(3.1) 


Using the transformation P — ^ hrst-order derivative term, induced 

by the background flow, can be eliminated. The results will be transformed back to the 
physical domain before making comparisons. The scattered near-held pressure is obtained 


by solving (3.1) using FEM and the far-held sound pressure is obtained by integrating 


the pressure distribution over surface, as described in (2.40) and (2.41). 


The commercial software COMSOL 4.4 is used to perform the FEM simulations. Sim¬ 
ulations are performed for a single serration and chord length c = Im, as shown in 
hgure For the boundary conditions, the normal velocity on the surface of the plate 
vanishes, while the pressure values are hxed over the wake half plane. In addition, as 
illustrated in hgure the walls on both the upper and lower sides of computational do¬ 
main represent Floquet periodic boundary conditions. The radiation boundary condition 
is implemented via Perfectly Matched Layers (PMLs), as shown in hgure The mesh 
is made of tetrahedron cells with quadratic shape functions. The mesh is highly non- 
uniform and is generated to accurately resolve the hydrodynamic pressure huctuations 
near the serrated edge and the acoustic pressure perturbation in the far-held. A mesh 
sensitivity test has been carried out to ensure the proper convergence of the simulation. 
In the hydrodynamic region, the mesh contains more than 10 grid points within one hy¬ 
drodynamic wavelength. The same ratio of 10 grid points per wavelength is used in the 
far-held, relative to the acoustic wavelength. The hnal mesh contains about 0.6 million 
elements. 

Figure 1^ shows the results for a straight trailing-edge at / = 1000 Hz and M = 0.1. 
The turbulent convection velocity is assumed to be Uc = 0.7U, where U = MqCq- The 
wavenumbers are ki = 2'Kf/Uc, ^2 = 0 and the amplitude of incident wall pressure 
gust Pi is unity. It can be seen that the PMLs do not cause spurious rehections. The 
solution was also compared with Amiet’s analytical solution. The maximum difference at 
/ = 1000 Hz was less than 1%. The error is partly due to the inhnite-chord assumption 
in the analytical model and in part to the FEM numerical errors. 
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-0.05 

H -0.1 

Figure 4. Computed pressure distribution scattered by a straight trailing-edge 


3.1.2. FEM model validation 

This section provides a comprehensive comparison between the analytical model for 


gust-induced far-held noise developed in Sec. 2.3.2 and the FEM model developed in 
Sec. 13.1.11 for different serrations and at different Mach numbers. The far-held sound 
pressure induced by a wall pressure gust of ^2 = 0 at diherent frequencies is chosen for 
comparison. Note that in the mathematical model, we use the second-order approxima¬ 


tion of the gust-response function C, as shown in (2.431. The chord length c is chosen 


to be unity with a span-to-chord ratio of 8. The Mach number and geometrical shape of 
the serrations can vary in different cases. 


Results are presented for a range of serrations, as illustrated in hgure 5(a) through 
|5(f)| In this study, we shall only focus on low Mach numbers, i.e. Mq ^ 0.2. The far-held 
pressure, obtained from the FEM model at 90° above the trailing-edge in the mid-span 
plane with Xz = 1, is plotted as 201ogiQ |p/(x,a;)| against the theoretical predictions. 
Results are provided for both the baseline (no serration, i.e. h = 0) and serrated trailing- 
edges. It can be found that a sound reduction of more than 20 dB can be achieved for 
the far-held sound induced by this specihc wall pressure gust. This reduction should 
however not to be confused with the sound reduction of the real trailing-edge noise, 
which comprises diherent wall pressure gusts at diherent values of ^ 2 . 

Figures [B(a)| to [5(^ show very good agreement between the theoretical and computa¬ 
tional results. The low frequency discrepancies originate from the inhnite chord assump¬ 
tion in the analytical model. The serration cases presented in hgures[5(a^to|5(d)|are not 


normally considered sharp enough, based on experimental observations ([Gruber 2012) 


to reduce the noise signihcantly. Figures [5(e)| and [5(f)| show the results at Mq = 0.2 for 
serrations with X/h = 0.5 and h/c = 0.1 and X/h = 0.3 and h/c = 0.1, respectively. It 
can be seen from these two hgures that for sharper serrations, the average error between 
the numerical calculations and theoretical predictions normally increases, which might 
be caused by the relatively slower convergence rate of the second-order approximations of 


(2.43) compared to that for the wide serrations. The agreement between the FEM results 


and the proposed model, however, is generally good, suggesting that the second-order 


solution does indeed give a reasonable good approximation for (2.43). The issue of the 
convergence of the iterative method will also be discussed later. 

3.2. The far-field sound spectrum 

A parametric study of far-field noise reduction was carried out by Howe (Howe|19916|a ), 


indicating the possibility of significant noise reduction, much higher than measured 


data ( 

Gruber 

2012 

). In this section, we shall use the second-order iterative model de- 

veloped in Sec 
Chase’s model 

2 and carry out a parametric study. For illustration purposes, we 
(Chase||1987 ) of the wavenumber spectral density. It is argued by 

adopt 

Ghase 


(1987) that the convection velocity Uc is weakly dependent on frequency and on aver- 
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(a) Mo = 0.1, \/h = 6, h/c = 0.025 (b) Mo = 0.2, \/h = 6, h/c = 0.025 




(c) Mo = 0.1, \/h = 3, h/c = 0.05 (d) Mo = 0.2, X/h = 3, /i/c = 0.05 




(f) Mo = 0.2, \/h = 0.3, /i/c = 0.1 


Figure 5. SPL at 90° above the trailing-edge in the mid-span plane with *3 = 1 due to a wall 

pressure gust of frequency ui with ^2 = 0 


age Uc « 0.7U. According to Chase’s model, the wavenumber spectral density is well 
approximated by 


n{uj,ki,k2) 


_ CmPlv*kl6^ _ 

{{ki - uj/U,)^6Uc/3v,)^ + {kl + kl)5-^ + x2)5/2 > 


(3.2) 
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where po is the density of air, and Cm ~ 0.1553, y « 1.33, r;* « 0.03C/. The turbulent 
boundary layer thickness 6 in (3.21 is approximated by ( Eckert fc Drake||1959 ) 

S/c = 0.382i^e-l/^ (3.3) 


where Rcc is the Reynolds number based on chord c. An inspection of (3.2) shows that 


the wavenumber spectrum peaks around ki = oj/Uc- Here, we assume ki = ui/Uc and 


obtain n(a;, ^ 2 ) by integrating (3.2) with respect to fci and then keeping the leading order 
terms ( Howe|[r9916 1, which yields 

-iCmPlvti^/UcYd^ 


n(a;, fc 2 ) 


(3.4) 


U,{{{u;/U,fAkl)5^+X^) 

Substituting (3.4) into (2.55) and using {pqv"IY[ d/C q) (Howe 19916) to non-dimensionalize 


the far field PSD yields 


where 'l'(x,a;) is defined as 


= —^(x w) 
{povlY{d/c^) 2 ^ ^ 


^ |/:(a;, 27 rm/A)|^ 




m——oo 


[{ujS/Uc)'^ + {211117:5/\Y + x^] 


(3.5) 


(3.6) 


The following figures are plotted using (3.6) for a variety of serration geometrical pa¬ 


rameters. Since much of the experimental work focused on the trailing-edge noise at low 
Mach numbers (Gruber et al. 2013), we shall only focus on the low Mach numbers, i.e. 


Mq ^ 0.2. The function £ in (3.6) is defined in (2.43) and we take the second-order ap¬ 


proximation here. Note that the incident pressure is also taken into consideration (Amiet 
|1978). The observer point is at 90° above the trailing-edge in the mid-span plane, namely 


(xi = 0, a ;2 = 0, X 3 = 1). It is worth pointing out that in figurej^both the far-held sound 
spectrum and the sound reduction spectrum are shown. 

The normalized sound power spectrum at Mq = 0.1 for different serrations are shown 
in hgures |6(a)| to 6(f)| The spectrum for a serrated trailing-edge with 4/i/A = 0.5 is 
shown in hgure |6(a) As expected, the sound reduction is approximately zero over the 
entire the frequency range of interest. Increasing the sharpness of the serrations gradually 
improves the sound reduction performance, as shown in hgures 6(b) 6(c )] and |6(d)[ For 
sufficiently sharp serrations, signihcant sound reduction is achieved over a wide range of 
frequencies, as shown in hgure [6(^ where the sharpness factor is 4/i/A = 10. The result 
obtained for a sawtooth serration with X/h = 0.2, h/c = 0.05 at Mq = 0.1 is shown in 
hgure |6(f)[ Comparing hgures |6(e)| and |6(f)| suggests that for already sharp serrations, 
further increasing the sharpness can provide a better high frequency noise reduction 
performance while the low frequency performance {kc < 10) remains unchanged. For the 
sharp serrations presented in hgure [6(f)] the far-held sound is reduced by about 10 dB at 
high frequencies. This better agrees with experiments where a noise reduction of up to 
7—10 dB is observed ( Dassen et a7|199^|Parchen et a7|1999 ). From hgure [6(a^ and [6(b) j 
it can be found that a slight noise increase may occur at low frequencies. In fact, the 
noise increase becomes even more pronounced at low frequencies when the Mach number 
is high, e.g. Mq = 0.4. The explanation of the noise increase at low frequencies will be 
given in Sec. [^ 

The results presented in hgure [^ were based on the second-order approximation. The 
convergence rate of different order solutions can be inspected by presenting the far- 
held sound spectrum using different-order approximations, as shown in hgure [^ where 
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Figure 6. The normalized spectrum of equation 3.6 (left axis) for straight and serrated trail- 
ing-edges, in solid and short dashed blue lines respectively, and the sound reduction spectrum 
ASPL (right axis) with Mq = 0.1 and an observer at 90° above the trailing-edge in the mid-span 
plane with xs = 1. 


the far-held spectrum using zero, hrst and second-order approximations are presented. 
Figure [^a) presents results for a wide serration with \/h = 2. As expected, due to the 
weak coupling between different modes the hrst and second-order solutions yield almost 
the same results. It is thus safe to assume that the second-order approximation gives 
an accurate solution for wide serrations. Figure [^b) shows the convergence results for 
a narrow serrations with X/h = 0.2. It can be seen that the difference between the hrst 
and second-order solutions is much smaller than that between the zero- and hrst-order 
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(a) A/fe = 2,/i/c = 0.025 


(b) \/h = 0.2, h/c = 0.05 


Figure 7. The normalized spectrum for straight and serrated trailing-edges obtained using 
different order approximations, Mq = 0.1, the observer is at 90° above the trailing-edge in the 
mid-span plane with *3 = 1 . 


ones. The maximum difference between the first and second-order approximations at high 
frequencies is less than 2 decibels. Thus, the second-order solution can be assumed to 
provide a reasonably accurate solution for narrow serrations even at high frequencies. 


3.3. Directivity patterns 

It is a well established fact that the trailing-edge noise directivity changes with fre¬ 
quency ( [Williams fc Hall 1970 Gruber et al. 20131. However, the effect of serrations 
on trailing-edge noise directivity has received very little research attention. Figures [^ 
and 1^ present the non-dimensional far-held PSD, (3.6), based on the second-order so¬ 
lution. Results are presented for straight and serrated trailing-edges with X/h = 0.4 at 
Mq = 0.1 and 0.4, respectively. As expected, using serrations has little effect on noise 
generation mechanism at very low frequencies, kc ^ 1. Results, however, show that the 
serrations can effectively reduce the noise at higher frequencies. As discussed earlier, this 
is believed to be primarily due to the destructive scattering interference effects. While 
most experimental investigations ( [Gruber et a/.||2M3 ) have focused on the capability of 
serrations for reducing the noise at small angles and 90 degrees above the trailing-edge, 
results in hgures [^ and [^ clearly show that serrations are more effective in reducing the 
noise at large radiation angles, i.e. towards the leading-edge, 9 > 90°. This is a very 
interesting result as noise measurement in laboratory environment is often limited to 30 
to 120 degrees due to the anechoic chamber room size constraint or reflection by the 
contraction nozzle, etc ( Gruber[[2012 Moreau fc Doolan[[2013 ). 

Results have also shown that the use of serrations can lead to significant changes to the 
directivity pattern of the scattered pressure field at high frequencies. While one would 
expect a cardioid pattern for straight edges at high frequencies, associated with the edge 
scattering of a half-plane ( Williams fc Hall|1970 |, results for serrated trailing-edges show 
that the directivity pattern is more dipolar with a clear peak at a specific angle which 
depends on both the serration sharpness and the Mach number. Numerical study of the 
directivity pattern for different serrations has shown that the expected cardioid shape 
gradually changes to a more dipolar shape as the serration sharpness increases and the 
directivity peak also gradually moves downstream, towards the trailing-edge. The dipolar 
behaviour of the noise from serrated trailing-edges means that the leading-edge region, 
9 = 180°, is much quieter than that for straight trailing-edge. Increasing the Mach 






































Prediction of noise from serrated trailing-edges 


21 


90 50 



90 60 


(a) fee = 1 
90 60 




(b) kc = 3 
90 60 


v30 


0 



(d) kc = 10 


90 50 



90 50 



Figure 8. The directivity patterns plotted against 8 in the mid-span plane {x 2 = 0 and 
r = \/ (xf x^) = 1) at Mo = 0.1 for serrations with X/h — 0.4 and h/c = 0.05. The far- 
held sound level in the hgures are lOlogj^Q • 10“^°). 
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90 50 



90 60 



90 60 90 60 




Figure 9. The directivity patterns plotted against 0 in the mid-span plane (x 2 = 0 and 
= \/ (a^i -|- x\) = 1) at Mo = 0.4 for serrations with X/h — 0.4 and h/c = 0.05. The far- 
held sound level in the hgures are lOlogj^Q (4'(x,ti;)/4 • 10“^°). 
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number also appears to move the peak angle towards the trailing-edge, as can be seen 
by comparing figures and It is also worth mentioning that in the case of high Mach 
numbers, see figure]^ the use of serrations can lead to considerable noise increase in the 
trailing-edge region (0° < 9 < 90°), for intermediate frequencies, 1 < fcc < 10. 


4. Comparison with Howe’s model 

The mathematical model and serration geometrical criteria developed by Howe have 
long been used as a tool to evaluate the effectiveness of trailing-edge serrations and 


estimate the level of noise reduction (Gruber 2012 Azarpeyvand et al. 2013 Jones & 
Sandber^|2012 1. However, it has repeatedly been shown that Howe’s model overpredicts 


the level of noise reduction (Dassen et al. 1996 Parchen et al. |1999[ Gruber 2012). To 
simplify the model, Howe assumes that the Mach number is sufficiently low to neglect the 
convection effect, the frequency are sufficiently high, satisfying ujh/Uc 1, the statistical 
property of the turbulence inside the boundary layer remains the same before and after 
passing the trailing-edge and the diffraction model is based on the Green’s function for 
straight trailing-edges and the slender-wing approximations. 

With the introduction of Ghase’s surface pressure wavenumber spectral density Model, 
Howe (]Howe|[r991a ) shows that the far-held PSD is given by 


Spp{u},x) 


iPovind/co){d/\^\y 


= C'm/7rsin^(6>/2) sin(^)4'/(a;). 


(4.1) 


where 


4'/(w) = 8{h/6y 


OO 
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(4.2) 

Even though the assumption of frozen turbulence is used in both models, Howe’s model 
differs from the model presented in this paper in several ways. In Howe’s model, the far- 
held sound pressure is based on a compact Green’s function. The Green’s function is 
obtained by making use of the slender wing approximation. The model developed in 
this paper, however, gives the scattered sound by solving the convected wave equation. 
Howe’s model neglects the effects of convection, so it is only valid at low Mach numbers. 
The new model is valid for any subsonic Mach number, as the convection effects have 
been incorporated in the convected wave equation. In addition, Howe’s model requires 
the high frequency condition uih/U 3> 1, while the new model requires kc > 1, as an 
inhnite chord is assumed in the derivation. 


Figures in this section represent the results obtained using Howe’s model, i.e. (4.1), 
and the model developed in this paper, i.e. (3.5). The normalized spectrum 'l'(x,a;) is 


defined in the same way as in (3.5). The result of the new model with the correction 
applied to the third term in (|2.46|), according to Amiet (1978), systematically increases 


the noise level by 6 dB. Figure 10'a) shows the noise prediction results for a serrated 


trailing-edge with \/h = 0.4, h/c = 0.05 at Mg = 0.1 using both models. The comparison 
of the results shows a clear difference between the two methods. At high frequencies, e.g. 
kc « 50, Howe’s model gives a sound reduction of about 13 dB while the new model 
predicts about 7 dB of noise reduction. Figure [T^b) presents the comparison for shaper 
serrations, with X/h = 0.2, h/c = 0.05 at Mq = 0.1. At kc « 50, the noise reduction 
predicted by Howe’s model and the new model are, respectively, 18 dB and 10 dB. It can 
be interpreted from the results that the new model provides a much more realistic noise 
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(a) (b) 



kc kc 


Figure 10. The normalized spectrum of Howe’s model and the new model, the observer is at 90° 
above the trailing-edge in the mid-span plane with xz = 1. (a) \/h = 0.4, h/c = 0.05, Mq = 0.1; 
(b) \/h = 0.2, h/c = 0.05, Mq = 0.1. 


reduction estimate and is more consistent with experimental observations 
1996 Parchen et a/.|[l9M Gruber|[2012 Gruber et aL||2013 |. 


)assen et al. 


It is very interesting to note that the zero-order solution accurately follows Howe’s 
solution at high frequencies, as shown in figure At intermediate frequencies, i.e. 2 < 
kc < 20, the zero-order solution oscillates strongly, but the mean value seems to be 
following Howe’s result. This is actually not hard to understood, as the Green’s function 
used in Howe’s model ( Howe|[r991a ), is in fact only valid locally. In other words, it does 
not include the coupling effect between adjacent sawtooth edges. Thus, at high frequency 
we expect the zero-order (without the coupling effect induced by the singular root and 
tip points) solution coincides with Howe’s results. The high-order solution, however adds 
the coupled interactions between different modes, and this coupling effect clearly reduces 
the sound reduction predicted by the zero-order solution at high frequencies. Thus, the 
large overprediction of Howe’s model is likely to have been caused by the choice of the 
Green’s function since the Green’s function is not able to take into account the coupling 
effects. 


5. Noise reduction mechanism 

In order to better understand the noise reduction mechanism, the pressure distribution 
over the flat plate surface is presented, see figures [m to As mentioned earlier, the 
incident pressure only raises the far-held sound by 6 dB systematically, thus it suffices to 
consider the scattered pressure distribution only. The scattered pressure, as mentioned 
in Sec. is essentially the pressure jump across the hat plate. As mentioned in Sec. |2.4[ 
the two non-dimensional parameters kih and kihe play an important role for effective 
sound reduction using serrated trailing-edges. In what follows, the scattered pressure 
distribution will be presented by hxing one parameter and varying the other. Note that 
the scattered pressure mentioned here is due to wall pressure gusts with /c 2 = 0. The 
discussion, however, also applies to gusts with k 2 ^ 0, as the streamwise number ki has 
the same value for different gusts. 

The scattered pressure on the hat plate is presented in hgure for different values 


of kih. The scattered surface pressure is obtained by evaluating the real part of (2.391 


using the second-order approximation (P, = 1) and then normalizing to unity. The results 
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Figure 11. The scattered surface pressure distribution at a fixed frequency for the same 
fci/ie = 7. (a) kih = 2; (b) kih = 4; (c) kih — 10; (d) kih = 20. 
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Figure 12. Scattered pressure on the serrated edge, (a) real part; (b) imaginary part. 
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presented in figure [TT] are obtained for kihe = 7, while kih varies between 2 and 20. The 
horizontal coordinate y' is normalized by the spanwise correlation length lyi and the 
vertical coordinate x' is normalized by the hydrodynamic wavelength Ai = 2'x/ki. Thus, 
the distance between the two adjacent vertical dashed lines corresponds to the spanwise 
correlation length, lyi. Figure [TT|a) plots the scattered surface pressure distribution near 
the trailing-edge for fci/i = 2. It shows that the scattered surface pressure field between 
two adjacent vertical lines is essentially in phase, so no strong phase variation within 
ly! occurs. Figure l^b) shows the scattered surface pressure distribution for kih = 4, 
and it can seen that little phase differences appear within a spanwise correlation length. 
Figure[TI|^ c) shows the scattered surface pressure distribution for kih = 10. It is clear that 
even if the spanwise correlation length becomes smaller, a pronounced phase difference 
still appears within adjacent vertical lines. Further increasing the value of kih to 20, 
as shown in figure |ll[ d) , decreases the spanwise correlation length, but enough phase 
difference still appears within the increasingly narrow ranges. 

To make the phase variation induced by the presence of serrations even clearer, the 
scattered surface pressure along the trailing-edge is presented in figure Each line cor¬ 
responds to a different value of kih. The real and imaginary parts of the pressure are 
shown in figure [T^a) and[T^b), respectively. The two figures are thus showing the pres¬ 
sure distributions at different instants. The red curve in figure [T^ a), which corresponds 
to the real part for kih = 2, remains almost entirely negative. The corresponding imagi¬ 
nary part, shown in red in figure [I^ b), has a phase which slightly changes signs over lyi. 
Since the signal oscillates between the real and imaginary parts, the phase chagnes sign 
only over a small fraction of the cycle. The black curves corresponding to kih = 20, on 
the other hand, show a strong variation within a spanwise correlation length in both fig¬ 
ure [T^a) and[I^b), indicating a strong phase variation over the whole cycle. Therefore, 
the phase differences of the scattered pressure are more likely to be strong and permanent 
for high values of kih. 

The scattered pressure distributions for different values of kihe are presented in fig- 
13 The values of kih is fixed at 10 while kih^. increases from 1 to 10 (see figures [T^a) 
d)). The horizontal and vertical coordinates are also normalized by lyi and Ai, re- 


ure 

to 


13 


spectively. Figure [I^a) presents the scattered surface pressure distribution for kihe = 1. 
As kihe is small, the distance between two adjacent dashed lines is very small compared 
to the serration wavelength. Thus, even though pronounced phase differences appear 
along the edge, one can hardly see any phase variations within a spanwise correlation 
length. The pressure distribution for kihe = 2 is shown in figure and no signif¬ 

icant phase variations are achieved. However, for fci/ie = 4 a clear phase variation of 
scattered pressure begins to appear within adjacent lines, as shown in figure [T^c). From 
figure [T^d), where kih^ = 10, it can be seen that pronounced phase differences appear 
within a correlation length lyi in the spanwise direction. To better visualize the phase 
variation appearing in the spanwise direction due to serrations, the scattered pressure 
along the serrated edge is shown in figurefor the same values of kihe as in figure [I^ 
Both the real and imaginary parts of the scattered pressure are presented. The tendency 
for large values of kihe to favor strong phase variations is clearly demonstrated. 

Finally, one is in a position to discuss the noise reduction mechanism by investigating 
the physical implications of the two parameters kih and kihe. From figure [TT] it is obvious 
that the criterion kih ^ 1 ensures an effective phase variation appearing along and 
near the trailing-edge in the spanwise direction. In addition, as shown from figure |13[ 
the condition kihe 3> 1 ensures that the phase difference appears within one spanwise 
correlation length, i.e. between two adjacent dashed lines shown in figure E] and |T^ 
Since the surface pressure within the regions bordered by dashed lines is coherent, and 
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Figure 13. The scattered surface pressure distribution for different serrations at a fixed 
frequency and kih = 10 (a) kihe = 1; (b) kihe = 3; (c) kihe = 6; (d) kih^ = 10. 
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Figure 14. Scattered pressure on the serrated edge for different serration geometries, (a) real 

part; (b) imaginary part. 
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since the range over which the phase difference appears is clearly much smaller than the 
acoustic wavelength, the far-held sound will be reduced due to destructive interference. 
Physically, this means that the phase differences induced on the hat plate in the spanwise 
direction due to the presence of serrations should be well situated within a correlated 


turbulent structure. Therefore, as demonstrated in both hgures and 1_3 the sound 
reduction is caused by the destructive interference of the scattered surface pressure due 
to the presence of serrations. 

Results in Sec. |3.2| have shown that in the case of wide serrations, a noise increase at 
low frequencies is also possible, especially at high Mach numbers. The reason is that at 
low frequencies, only little phase variation is induced by the presence of serrated trailing- 
edges in the spanwise direction, but the wetted length of the trailing-edge is in fact 
much longer than that of a straight one. Thus, the net effect of phase interference can 
be constructive, which leads to an increase in the far-held noise. This is more likely to 
occur when kih is small, i.e. Mq is large, frequency is low or the serration is wide (small 
value of h). 

Based on the preceding discussions, one can conclude that the minimum effective 
serration length required for noise reduction can be obtained from hmin = min(h, he). It 
is then straightforward to combine the two conditions stated above into one, i.e. kihmin ^ 
1. Therefore, in order to achieve an effective noise reduction in the far-hied, the geometry 
of the serrations should satisfy kihmin ^ 1- Based on the interference results in hgurefTT] 
and 13 it can be found that a common rule of thumb is kihmin ^ with higher value 


more favourable. 


6. Conclusion 

A new mathematical model is developed in this paper to predict the sound radiated 
by serrated trailing-edges. The model begins with establishing an idealized scattering 
problem, resulting in a mixed boundary value convective wave equation problem with 
complex boundaries. This leads to a set of coupled partial differential equations, which 
cannot be solved using the standard separation method. A solution is obtained based on 
Fourier expansion to separate the variables and Schwarzschild’s method together with 
an iterative technique to solve the resulting coupled equations. The far-held sound is 
evaluated using the surface pressure integrals. The PSD of far-held sound is related to 
the wavenumber spectral density of the wall pressure beneath the turbulent boundary 
layer near the trailing-edge using Amiet’s approach ( Amiet|[l976fc 1978). 

The results obtained using the new model agree well with FEM computations, sug¬ 
gesting that the model developed in this paper captures the scattering process and gives 
correct predictions for the sound generated by serrated trailing-edges. It is shown that 
the coupling effect must not be ignored and as a result the new model can predict the 
sound reduction more accurately than Howe’s model (Howe 19916|a |. The results ob¬ 
tained using the new model agree better with experiments, in which the average sound 
reduction is reported to be up to around 7 dB. The directivity results show that serra¬ 
tions can significantly reduce the noise in the area near the leading edge and that at high 
Mach numbers the use of serrations can lead to noise increase at small angles. 

The physical mechanism for noise reduction is found to be interference effects in the 
wall pressure fluctuations due to the presence of serrations. Two non-dimensional pa¬ 
rameters are found to be critical. First, kih ^ 1 to ensure the existence of strong phase 
variation in the spanwise direction. Second, kihe^ I to ensure that the phase differences 
along the edges are correlated in the spanwise direction. The sound reduction generally 
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increases as the serration sharpness increases, but if the serrations are already sharp 
enough, further increasing the slope only affects high frequencies. 

The results obtained using Chase’s turbulent boundary layer spectrum model do not 
appear to explain the noise increase observed in experiments at high frequencies ([Parchen 


et al. 1999 Oerlemans et al. 2009 Gruber 20121. This suggests that the wavenumber- 


frequency spectra of the surface pressure fluctuations (see Sec. 2.41 are not accurate or 
that this noise increase is due to some other mechanisms, such as the high intensity flow 
through serration valleys. In addition, the current model assumes perfect correlations 
in the streamwise direction, which may not be sufficiently accurate according to the 
measurements of Gruber (2012|. Thus, the current model may be further improved by 
incorporating more physical parameters such as the streamwise correlation length. 
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Appendix A. 

A.l. Scattered pressure of second iteration 
The solution after the second iteration can be expressed as, 

0) = Nix) + (a;) + (x), (A 1) 

where Nix) and C^^^x) are defined in Sec. and C^'^\x) whose entry corresponding 
to mode n' is 
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A.2. Far-field sound pressure of second iteration 
The function involved in the second iteration can be defined as: 
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The second iterated solution falls into the same pattern, 
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